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I .  INTRODUCTION 

The  thermoelasticity  problems  found  in  mechanics  usually  revolve 
around  one  of  four  areas : 

1.  The  measurement  of  a  body's  surface  temperature  from  its 
surface  deformation  field. 

2.  The  measurement  of  thermal  gradients  at  the  surface  of  a  body 
from  its  surface  deformation  field. 

3.  The  measurement  of  thermal-induced  deformation  at  the  surface 

of  a  body. 

4.  The  measurement  of  thermal-induced  stress/strain  at  the 
surface  of  a  body. 

The  objective  of  this  report  is  to  present  the  theory  and  experimental 
laser  speckle  interferometry  techniques  that  were  used  to  obtain  these 
measurements.  Section  II  lays  the  foundation  j.-r  using  the  laser  to 
make  these  noncontact  measurements.  Since  laser  speckle  interferometry 
is  used  to  measure  in-plane  motion  only,  certain  restrictions  have  to 
be  imposed  on  the  theory. 

Although  the  theory  is  developed  for  three  dimensions,  necessary 
simplifications  are  used  to  reduce  this  theory  by  one  dimension. 

Section  III  consists  of  two  sample  problems  that  were  used  to  test  the 
theory  of  Section  II.  The  thin  circular  flat  plate  and  uniformly 
heated  rod  are  treated  in  Section  II.  Section  IV  presents  the 
conclusions  drawn  from  this  work. 

II.  THEORETICAL  ANALYSIS 

A.  General  Thermoelasticity  Theory 

The  linear  theory  of  elasticity  may  be  used  to  predict  the 
deformation  field  of  a  body  resulting  from  thermal-induced  stress  and 
strain.  A  temperature  change  at  an  arbitrary  location  in  a  body  may  be 
predicted  from  the  local  deformation  field.  The  theory  to  follow 
assumes  that  the  body  consists  of  an  isotropic  hookean  material  which 
obeys  linear  thermoelastic  behavior.  For  anisotropic  materials,  the 
theory  becomes  significantly  more  involved  and  will  not  be  treated  in 
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this  analysis.  Much  of  the  analysis  will  be  presented  in  tensor  form 
to  simplify  the  work. 

The  strain  displacement  relations  from  the  theory  of  elasticity  may 
he  expressed  as 


e 


ij 


[3Ui 

*!24- 

3U . 

6  — - 

l8*J 

8*i| 

ij  3Xj 

(i) 


Equation  (l)  expresses  the  strain  field  to  the  displacement  field 

U  in  the  rectangular  coordinate  system  The  equilibrium  equations, 

i  J 

which  relate  the  stress  field  a.  .  and  internal  body  forces  ,  may  be 

^  J 

expressed  as 


3a-  • 

_^B.  =  0 


(2) 


The  following  analysis  was  obtained  from  reference  [1].  The  boundary 
condition  equation,  which  relates  the  stress  on  an  arbitrary  plane 
in  a  body  to  the  body's  internal  stress  field,  is  given  as 
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where  n  are  the  direction  cosines  of  the  plane.  In  expanded  form,  the 
i 

components  of  Eq.  (3)  become 
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For  a  body  free  of  thermal- induced  strain,  the  stress/strain  relations 
may  be  expressed  as 


exx  ~  E  txx  ~  y^ayy  +  azz3 

eyy  =  E  [°yy  ' 

czz  =  E  [°Zs  '  *  “yy^ 

Exy  =  G  °xy 

eyz  “  (T  °yz 
1 

ezx  G  °zx 


(4) 


In  Eq.  (4),  E  is  the  modulus  of  elasticity,  u  is  poissons  ratio,  and  G 
is  a  Lame'  constant. 

Using  the  general  stress/strain  relations  and  the  concept  of  thermal- 
induced  strains,  the  change  in  temperature  AT  in  some  local  region  G  of 
a  body  may  be  related  to  the  displacement  field  in  the  region  G  [1]. 
First,  the  thermal-induced  strains  for  a  body  uniformly  heated  and 
subject  to  a  free  expansion  is 
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(5) 


et  =  aAT 


where  is  the  local  normal  strain  in  any  direction  t  of  ^  a  is  the 
coefficient  of  thermal  expansion,  and  AT  is  the  temperature  change  over 
ft.  The  temperature  coupled  stress/strain  relations  are  [1] 
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If  all  the  stresses  are  temperature  induced  for  a  differential 
element  with  only  thermal -induced  strain,  then 


EaAT  =  a  -  p  (a  +  a  ) 
x  v  y  z 7 


EaAT  =  a  -  p(o  +  a  ) 
y  v  x  z' 


EaAT  =  a  y(a  +  a  )  ; 

z  v  x  y'  * 


(7a) 

(7b) 

(7c) 


then 


3EaAT  -  (ox  +  ay  +  az)(l  -  2y) 


(8) 
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Subtract  Eq,  (7b)  from  Eq.  (7a)  and  Eq.  (7c)  from  Eq.  (7b)  to  obtain 


°x  "  ~  °y  +  =  0  ^  °x  =  CTy 

°y  -  V°z  -  °z  +  vay  =  0  “*  °y  =  °z 


or 


a  =  a  =  a 
x  y  z 


EaAT 

1  - 


(9) 


Therefore,  expansion  of  the  differential  element  may  be  prevented  by 
the  application  of  a  hydrostatic  stress: 


-EaAT 

at  1  -  2p 


(10) 


The  equilibrium  equations,  Eq.  (2),  without  body  forces  present  take 
the  form 
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Consider  a  differential  stress  element  with  the 
given  in  Figure  1,  The  stress  vectors  on  the  faces 


unit  normals ,  as 

n  ,  n  ,  and  n  are 
x*  y  z 
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Figure  1 


Differential  stress  element  with  unit  normals  nx>ny>nz* 
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(12) 


S  =  a  i  +  a  j  +  a  k 
x  xx  xy  u  xz 


S  =  a  i  +  a  j  +  a  k 
y  xy  yy  yz 


Sz  =  “xs  1  *  ays  J  *  °zz  k 


respectively.  The  equilibrium  equations  may  now  be  expressed  as 
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(13) 


The  stress  components  may  now  be  written  in  terms  of  strain  components . 
Let 
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In  Eq.  (15),  t  is  the  deformation  vector  between  the  loaded  and 
unloaded  conditions  for  the  solid  and  is  given  by 


=  ui  +  vj  +  wk 


(16) 


Substituting  Eq.  (15)  into  Eq,  (12)  yields 


+  2s. )  +  i 


sx  =  G  VVU  +  &  I  +  1  v  ‘  p 


[  ->■  9d  1  .  A.  2yG  ^  ^ 
sy  =  G  lvv  +  W  +  J  v  p 


(17) 


Now,  with  thermal  stress  persent,  a  -l  term  must  be  included  for 

each  S.  stress  vector.  Therefore, 
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Now,  for  equilibrium 
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or  V  •  S 


V  •  S 


V  •  S  =  G 


V2v  + 


V2W  + 


1  9  n  . 

Ea  9AT 

1  -  2p  Bx  V  p 

1  -  2y  3x 

1  9  v  .  D 

Ea  9AT 

1  -  2y  By  P  j 

1  -  2y  9y 

1  9  =*  -*•  1 

\  Ea  9 AT 

I  -  2y  9z  V  P 

|i  1  -  2y  3z 

=  0 


=  0 
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rearranging  terms  yields 


2  1  I  -»  En  -*• 
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(21) 


let 


X  = 


G  +  'yKi -"M 


(22) 
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The  deformation  vector  function  p  may  now  be  expressed  in  terms  of  the 
strain  potential  as 


2Gp  = 


V<j> 


(23) 


Since  VV 


VV 


7  v,  then  Eq.  (2l)  may  be  written  as 


G 


V<)> 

2G 


2-x 

V  V  + 


1 

1  -  2y 


£ 

2 


2  -  2y 

=  TTT2I7 


VV 


2  $ 
7 


(24) 
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a  sufficient  condition  is  that 
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An  equivalent  expression  for  Eq.  (24)  is  [2] 


(A  +  °)  ||  +  GVU  -  (3A  +  2G)a  =  0 


(A  +  G)  ||  +  Gv2v  -  (3A  +  2G)a  -  0 


(A  +  G)  ||  +  Gvw  -  (3A  +  2G)a  ||i  =  0 


(29) 
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where  e  is  the  dilatation,  or 


9u  +  9v  +  8w 
9x  9y  9z 


(30) 


The  relationship  between  the  stress  field  and  temperature  change  may  be 
expressed  as  [2] 
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B.  Laser  Speckle  Interferometry  Theory 

Laser  speckle  interferograms  are  most  commonly  used  to  make 
surface  displacement  measurements  of  a  deformed  body,  Figure  2 
illustrates  the  basic  method  for  making  a  laser  speckle  inter ferogram. 
When  the  diffuse  surface  of  a  structure  is  illuminated  with  coherent 
radiation,  a  grainy  speckle  effect  is  imaged  by  the  eye  or  film  plane 
of  a  camera  due  to  the  interference  of  light  from  the  structure.  This 
speckle  effect  is  enhanced  when  the  structure  has  microscopic  surface 
irregularities.  If  the  optical  configuration  remains  fixed,  the 
speckle  pattern  of  the  test  object  may  be  recorded  on  the  film  plane  of 
a  camera.  Further,  if  the  structure  is  deformed,  the  speckle  points 
shift  with  the  deformation  and  a  second  exposure  of  the  deformed 
speckle  pattern  can  be  made. 

Using  a  technique  of  double  exposure,  spec^.  e  interferograms  of  a 
structure  are  normally  made  by  photographing  the  speckle  pattern  in  a 
reference  and  deformed  configuration.  A  beam  of  laser  light  is  then 
passed  through  a  region  of  the  double  exposure  where  the  local  deforma¬ 
tion  is  desired.  As  the  beam  passes  through  the  film,  the  deformed  and 
undeformed  speckle  recorded  there  diffract  the  laser  light  and  cause  an 
interference  effect  on  a  viewing  screen.  A  diffraction  halo  modulated 
by  light  and  dark  bars  of  light  is  produced  where  the  distance  2d 
between  bars  is  inversely  proportional  to  the  distance  between  the 
undeformed  and  deformed  speckle  on  the  film  plane.  A  normal  to  the 
light  and  dark  pattern  indicates  the  axis  of  deformation  of  the  speckle. 
The  theory  to  be  presented  assumes  that  the  deformation  region  illumi¬ 
nated  by  the  laser  beam  in  reconstruction  is  uniform  and  that  the  linear 
optical  theory  is  applicable. 

Figure  3  illustrates  the  reconstructed  diffraction  halo  modulated  by 
light  and  dark  bars  of  light.  From  the  linear  theory  [3],  the  displace¬ 
ment  in  the  0  direction  of  a  point  on  the  body  is  given  as 


SXf 

■0  2d 


(32) 
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(A)  FORMATION  PROCESS 


(B)  RECONSTRUCTION  PROCESS 


Figure  2.  Laser  speckle  interferometry  configuration. 
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Figure  3.  Diffraction  halo  geometry. 
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where 


S 

X 

f 

d 


film  scale  factor  (magnification  ratio) 
wavelength  of  laser  illumination  source 
distance  from  interferogram  to  analyzer  screen 
distance  from  central  bright  spot  to  first  minima 
displacement  of  the  point  illuminated  by  the  laser  on  the 
object  in  the  0  direction 


The  vertical,  U  ,  and  horizontal,  Ujj,  components  of  displacement 
may  be  obtained  from  U0  as 


UR  =  U0  cos  0 


S\f 

2d" 


cos  0 


U  =  UQ  sin  0 
v  0 


Sxf 

2d" 


sin  0 


(33) 


and  from  the  geometry 


d 


COS  0 


d 


T 

V 


sin  0 


Therefore, 


(34) 


U, 


H 


SXf 


u 

v 


sxf 

2d~ 


(35) 


Laser  speckle  interferometry  can  be  used  to  make  very  accurate 
measurements  of  the  in-plane  deformation  of  solids.  However,  out-of 
place  deformation  cannot  be  usefully  measured  with  this  technique. 
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Therefore  the  technique  suffers  slightly.  In  many  cases  (i.e.,  thin 
plates),  out-of-plane  deformation  estimations  can  be  made,  which  allow 
accurate  measurements  of  temperature  change  to  be  predicted.  In  a 
typical  application,  laser  speckle  interferometry  may  be  used  to 
estimate  the  u  and  v  components  of  deformation,  a  postulation  between 
the  dependence  between  the  w  component  of  deformation  and  temperature 
change  may  be  made,  and  the  temperature  change  of  the  solid  may  be 
estimated  based  on  the  simplification  of  Eq.  (29). 

C.  Thermal  Expansion  of  a  Heated  Rod 

For  a  uniformly  heated  rod  (Figure  4),  the  expansion  with 
temperature  change  may  be  estimated  from 


Ae  =  aLAT 


(36) 


where  Ae 
a 
L 
AT 


change  in  the  rod  length 
coefficient  of  thermal  expansion 
rod  length 

temperature  change  of  rod  (uniform) 


For  a  12-inch  rod  with  a  coefficient  of  thermal  expansion  equal  to 
1.66  x  10-5  °CT1,  a  .502  °C  temperature  change  is  required  to  achieve  a 
Ae  =  .0001  inch.  Figure  5  illustrates  the  dependence  between  Ae  and  AT 
for  various  L. 


Figure  4.  Thermal  expansion  of  a  uniformly  heated  rod. 
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0  .001  .002  ,003  .004  ,005 

Ae  ( in. ) 


-5  -1 

Figure  5,  AT  versus  Ae  for  various  rod  lengths  (L),  a  =1.66  x  10  °C  . 
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D.  Thermal  Stress  Deformation  of  a  Thin  Circular  Plate  Neglecting 
the  z  Component  of  Deformation 

The  general  temperature/deformation  equations  take  the  form 


Ea 

1  -  v 


AT 


(37) 


The  Laplacian  and  Gradient  in  cylindrical  coordinates  are 


2 

V  <|> 


=  1 2_  (r  2iU  i_  ift  *  it 

r  8r  1  8r)  777  77 


5*  =  21  "e  *  I  21  e  *  |±  ; 
v<p  gr  r  r  ge  e  az  z 


(38) 


For  a  circular  plate  with  a  radial  temperature  distribution  AT  =  AT(r) 
and  neglecting  the  z  component  of  plate  deformation, 


1  3_ 
r  3r 


r|* 

ar 


Ea 


AT 


<f>  = 


I  -  y 

<j>(r)  only. 


(39) 


For  such  a  symmetric  circular  plate. 


V(j> 


ar 


e 

r 


(40) 
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Therefore 


1  3 


111  = 

Ea  AT 

3r  j 

l  -  y 

111  = 

rEa  AT 

3r  1 

1  -  V 

r  11  = 

fT  rEa  , 

r  3r 

1 

( — 1 

9<J>  _ 

9r 

1  fT  rEa 

r  J  1  - 

A 

AT  dr  +  C, 


cr 

AT  dr  +  — 


C1  =  0  since  p  =  0  at  r  =  0 


.4.13*:- 


7G  ~  2G  3r 


(41) 


and 


[  1  f  r  rEa 

s  \zgzj 0  r^nr 


AT  dr  e 


(42) 


The  radial  deformation  measured  with  laser  speckle  interferometry  may 
he  expressed  as 


mAfS  -*■ 
pr  er  ~  6r 


(43) 
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where  S  =  film  scale  factor 

f  =  viewing  screen  to  interferogram  distance 

X  =  light  source  wavelength 

m  =  fringe  order 

x(r)  =  fringe  spacing 

Substituting  Eq.  (43)  into  Eq.  (42), 


roXfS  A 

xCr7  ep  ^  jo  i  -  M 


1  f  1r^-  sTdr  e_ 


or 


=  (sects?  I  rrv mr} 


x(r)  = 


Is  (A  p)I 


mXfSr 


( iMjf 


rATdr 


x(r) 


.  M 


mXfSr 


/r 

rAT(r)di 


(44) 


(45) 


For  laboratory  analysis,  x(r)  is  measured  as  the  first  minima  to  first 
minima  spacing;  i.e.,  m  =  1.  Therefore, 


(46) 
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E.  Thermal  Stress  Deformation  of  a  Thin  Circular  Plate 
Correction  for  the  z  Component  cf  Deformation 

Consider  a  thin  circular  plate  as  shown  in  Figure  6 
potential  function  for  this  case  is  ij>  =  f(r,z).; 
therefore 


1 

2$ 


+ 


9<f>  -»■ 
TZ  ez 


> 


and 


1  9<j> 

r  Sr 

N 

1  84> 

r  "Sr 

9jfr 

9z2 


a2<t> 

77 


Ea 

1  -  y 


AT 


To  simplify,  let 
AT  =  AT(r); 

thus 


13K  92<t>  .  92<t> 

7T  7? 

9r  9z 


Ea 

T-TJ 


AT(r) 


and 


P  = 


|9<f>  A  ^  9*  2  I 
TSS  =  25  11?  6r  9l  6z| 


V(J>  _  1 


mXfS  '  9<t> 

xC.r)  r  2G9z  z 


with  a 

.  The 

(47) 


(48) 


(49) 

(50) 
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Figure  6. 


Thermal  stress  deformation  of  a  thin  circular  plate 

with  p  and  p  components  of  deformation, 
r  z 
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or 


1  3<(>  _  mAfS 

7G  9r  xC^J 


it  =  2GnafS 

9r  x(r) 


(51) 


The  thermal  profile  equation  then  takes  the  form 


1  2GmAfS 
r  xCr) 


rrv  4TW 


(52) 


With  laser  speckle  interferometry,  the  e  component  of  displacement 

z  2  2. 

cannot  be  measured.  Only  an  estimate  of  the  significance  of  9  $/9z  in 
the  thermal  profile  equation  can  be  made.  Assume  that  the  plate  is  of 
thickness  2t,  which  is  thin  enough  that  there  is  no  temperature  variation 
in  the  z  direction.  Figure  7  illustrates  the  assumed  deformation  pro¬ 
file.  From  this  figure,  the  deformation  is  linear  with  z  and  takes  the 
form 


P 


z 


=  2GaAT  z 
9z 


(53) 
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plate  deformation 
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& 
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-p 
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<D 
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Ph 
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Ph 

a) 

d 
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d 

< 

a 

<D 

N 

(U 
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B 

hD 
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in  the  z  direction  of  a  heated  flat  plate. 


Substituting  Eq.  (53)  into  Eq.  (52)  yields 


1  2G  mXfS  +  2GmXfS 
r  X\r) 


h  (xw)  *  2GoAT(r)  ■  r=TT 1T(r) 


and 


W(r)  |jS2- 


-  2 Gal 


(54) 


Define 


«r)  =  i 


8r  WJ\ 


B  = 


mXfS  Jl  -  yl 
2a  1  y 


(55) 


then 


AT(r)  =  BD(r) 


(56) 


Let 


y(r) 


=  xtrj 


then 


AT(r) 


(57) 


(58) 
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F.  Thin  Plate  Thermostress  Equations  in  Rectangular  Coordinates 


In  rectangular  coordinates. 


op-*’  -  9<f>  9<f>  ~  9 <f)  * 

2Gp  =  V<f>  =  e  +7rLe  +  *1-  e 
9x  x  9y  y  9z  z 


V2<J>  -  AT 

1  ^  u 


2  2  2 


3x  3y  3z 


(59) 


By  substituting, 


9 


■tsx  <2Qpx) *  ly  (2V  *  5S  (*»>»)  -  rrr41 


3px  *  !i .  !!* 

9x  9y  9z 


=  a 


1  +  yl 

T-nti 


AT 


(60) 


In  general, 


1  , 

l  -  y 

9p  9p  9p  ^ 

_£  +  _V  +  A 

a 

1  +  y 

9x  9y  9z  I 

(61) 


For  thin  plates,  there  are  two  independent  approximations: 


1.  p  =  aATz 

Zi 


^2 

9z 


=  aAT 


9pz 

2.  -JL  8  0 

9z 
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It  is  important  to  note  that  in  many  instances,  Eq,  (6l)  will  not 
satisfy  all  the  boundary  conditions  of  a  problem.  For  approximation  1, 


9p  9p . 


-  -  r  Nr  ir  *  ^ 


Mnr" 


(62) 


Let 


-  I- 

g  \l 


-  V 


1  - 


+  V 

T^y 

r“nr 


(63) 


Therefore 


9PX  9Py 
AT  "  ^l  lax-  +  ay 


( 64) 


For  approximation  2,  let 


then 


1  1  -  y  1 

a  |1  +  y| 


(65) 


aT  _  .  Japx  +  apy 
aT  -  ^2  fer  +  aF 


(66) 


Using  laser  speckle  interferometry,  approximation  1  may  he  expressed  as 


AT 


lJ^SAf 

2 


and  approximation  2  is 


AT 


^2SAf  fa_  (l_ 

| ax  ydR 


(67) 


( 68) 


32 


G.  Numerical  Differentiation  of  the  Thermal  Induced  Deformation 
Field  to  Predict  Temperature  Change 

Suppose  the  temperature  change  upon  heating  or  cooling  is 

desired  at  some  location  i,j  on  a  body.  The  x  components  IL, 

U.  n  and  the  y  components  V.  , ,  V.,  V.  ,  of  displacement  have  been 
i+l  j-1  J  J+1 

determined ,  Figure  8  illustrates  the  geometry. 


y 


* 


T 


i 


v . 
J 


+1 


MEASUREMENT  POINT 


DISPLACEMENTS  ARE  OBTAINED 

FROM  THE  LASER  SPECKLE 

INTERFEROGRAM  AT  LOCATIONS 

i-1,  i,  i+l,  j,  J-l,  J+l  AND 

USED  TO  DETERMINE  AT... 

ij 


Figure  8,  Geometry  for  measuring  temperature  change 
AT  at  location  i,j. 
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The  approximation  to  the  derivative 


(  x  3p*  +  3py 

p(x,y)  =  —  +  ay- 


may  he  obtained  as  follows: 


3p 


,  U.  -  U.  ,  U.  ,  -  U 
x  „  1  |  i  i-l  x  1+1 _ 1 

x~  =  ? — a  d 


u.  ,  -  u.  , 

1+1  1-1 

2d 


9p  V.  -  V.  V.  .  -  V. 

y  «  i  \  j  tL  +  J+1,  J 
w  '  ? I a  a 


V.±1  -  V.  . 
J+1  J-l 

2H 


or 


p(x,y)  * 


U.  n  -  U.  ..  +  V.  .  -  V. 
1+1  1-1  J+l  J-l 

2d 


(69) 


(70) 


Using  laser  speckle  interferometry. 


p(x,y) 


SAf 

XT 


(71) 


H.  Deformation  of  Circular  Flat  Plates  in  Cylindrical 
Coordinates 

The  following  set  of  equations  are  derived  in  cylindrical 
coordinates  to  relate  the  change  of  temperature  field  in  thin  circular 
flat  plates  to  the  corresponding  deformation  field.  For  cylindrical 
coordinates, 
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G 


(72) 


.  „2  1  -»-»  l  V<J> 

17  *  rrm  V7-[  g 


Ea 


1  -  2P 


VAT 


^  1  *  1  9<J>  *  .  3<f>  A 

dr  v  v  1-kT  0  9z  z 


,,2  -  18 

7  *  '  F5r 


dip 

r  Sr 


2  2 
.  1  3  **  3  <(> 

~~2  — 2"  +  — ? 

r  30  3z 


V-*  -  FlFK^Flr^)*!^) 


The  derivatives  of  the  cylindrical  vectors  e0  and  er  are  [4] 


30  e0  er 


—  e  =g 
30  r  0 


Since 


->■  _  Vd> 

P  =  2G  ’ 


then 


12  1  -*■-* 
G  v  +  rr-^r  vv* 


p  - 


Ea 

rr^r 


}  VAT 


2  1  -*-* 
V  +  TTT-ZZ  vv* 


P  = 


(l  +  y)a 

TV  2y 


VAT 


(73) 

(74) 

(75) 

(76) 

(77) 

(78) 

(79) 

(80) 


35 


Let 


y  = 


l 

1  *-  2p 


(81) 


B  = 


2  (l  +  y) 

1  -  iy 


(82) 


Now  Eq.  (80)  may  be  written  as 


|V  +  yVV* 


p  =  BVAT 


(83) 


Expanding  Eq.  (83), 


1  3 

r  9r  |r  9r 


2  2 

+  in-  -w  +  —  n  lip  e  +  p„e.  +  p  e 
~~2  72/  \  r  r  ^6  0  yz  z 

r  90  9z 


+  Y\9r  6r  +  r  96  6  6  9z  6z)  \9r  6r  r  90  60  9z  z 


p  e  +  p„e„  +  p  e 
^r  r  1  0  z  z 


S  I  —  e  +  i.  —  e  +  —  e  )AT 
®  '  9r  er  r  90  e0  9z  z / 


.  (84) 
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Further  expansion  yields 


In  cylindrical  coordinates, 


pre0 


2 

3  *  8p  A 

^2”  6r  +  90”  e0  +  30“  60  ~  prer 


(86) 


(87) 


8p0  - 
90  e0 


p0er 


(88) 


2 

^  Pq  A  0  A  ^^0  ^  * 

g02  60  90  er  90  er  ~  p060 


(89) 
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Substituting  Eqs.  (86),  (87),  (88),  and  (89)  into  Eq.  (85)  yields 


1  a_ 

r  9r 


e 

z 


1  ^  pr  *  2  ^Pr  A 

+  “2  ~r  er  +  “2  ~  e0 
r  96  r 


+ 


/2 

i  3  pe 


r  90 


2  9pe  A  pe  A 
T  90"  er  '  T  e0' 

r  r 


+ 


e 


0 


+ 


e 

z 


+  yL  {li¬ 
ar  lr  3r 


l  a 

+  F3F  pe 


e 

r 


+  Ii_ 

r  90 


(li. 


,  l  a 

F3FP@ 


+ 


a_  ] 

az  pzj 


+ 


l  a 

r  F0  P0 


e 

z 


=  e 


(aAT  A  1  3AT  A 
~  er  +  r  W  e0 


(90) 
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Assume  for  the  cylindrical  plate  that  due  to  symmetry,  Pg  =  0. 
Therefore , 


Assume  that  =  p^(r)  and  p^  =  pz(r,z);  then 


1  9_ 
r  ar 


e 

r 


+  II- 

.  r  3r 


+ 


e 

z 


(92) 


From  Eq.  (92), 


3  AT 
36 


=  0 


(93) 


Examine  the  e  variation  in  Eq.  (92); 
r 


13  I  apr  I  ;  prer 
F37  r  —  er"~2-  +  Y 


I?  {?  lr  W  *  5^}  ‘er  =  6 


(94) 

ar  r 


or 


1  3  J  9prl  Pr 

7  37  f  — )  "  7  +  Y 


3  (l  3  (  \  ,  9pzl  _  0  3AT 

3r  |r  3r  \Pr/  3z  j  3r 


(95) 


Expanding  Eq.  (95), 


8pr  pr 


32pr 


lli  +  i!!£-!£+y  —  +  —)  +  y  — 

r  3r  T  '  3r  hr  r  '  3r3z 
3r  r  l  / 


=  6 


3AT 

3r 


(96) 


32p  3p  p 

r  +  x  r  Mr 


9\ 


2 

9  Pr  x  9Pr  Pr 

TT  +  7  3? - ?  +  Y  “T"  +  r  W-  -  Y  T  +  Y  3r3z  "  3r 

3r  r  3r  r 


=  8  2*1  (97) 


32P. 


3r  r 


Pr  -  T  lii  =  6  ^  (98) 

Y  3r3z  P  3r  x  ' 


(1  +  y) 


32p 


9r 


r  1  9pr  pr 

"  +  r  9r - ? 

r 


+ 


32p 

z 

3r3z 


=  3 


3AT 

3r 


(99) 
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In  general. 


assume  p  =  aZAT(r),  Then 


•2p. 

Y  3r3z 


and 


(1  +  y)1 


;  Y^CcxAT) 

=  1 

'  2 

8  p  n  9p 

r  +  1  Kr 

pr 

~  2  r  3r 

"  “2 

3r 

r 

.  3pr  °rl 

U 

3AT 

3r 


«  (3  ~  ay) 


d2pr 

Tf+  r  3r - 2 

3r  r 


3AT 

3r 


8  -  ay  I  3AT 
1 1  +  Y  I  3r 


From  Eq.  ( 102 ), 


B  ^  ay 

!  +  Y 


2a  +  2ya  a 
l  -  2'y"  "  T1T2Z 

- T - 


=  r 


-l 


1  + 


r^r 


r  =  l  rHp 


Equation  (102 )  may  now  "be  expressed  as : 


(100) 


(101) 


(102) 


(103) 


(104) 


(105) 


(106) 
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and 


iL  (il  ^  pr)  _  9AT 
9r  |9r  r  I  9r 


Finally, 


r 


9r  r 


+  c 


AT(r) 


(107) 


(108) 


Equation  (108)  is  based  on  the  following  assumptions: 


AT(r) 

Pr(r) 

0 

aZAT(r) 

I.  General  Thermostress  Equations  for  Thin  Plates  in  Rectangular 
Coordinates 

The  general  thermostress  equations  in  rectangular  coordinates 
have  the  form 


AT  = 
Pr  F 
pe  = 

pz 


G 


+ 


1 

I  V  3T 


Ea 

1  -  2U 


VAT 


or  if 


(109) 


e 


1 

tt-^t 


Ea  _  a2  (l  +  y) 
7  =  =  1  -Tu 


then 


(7 


2 


+ 


-A**  \ 

yW* ) 


BVAT 


(no) 
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Expanding  Eq.  (110)  yields 


2  2  2 

.3  ^  3  .  3 

lax  ay  az 


(a  A  I?)'*  ^  A 

pxi  +  Pyj  +  Pzk)  +  1  |a^  1  +  a7  j 


a  ;l  fpx  +  3py  .  9pz 
+  Sz  k  E"  +  3y  W 


0  J  3AT  t  .  3AT  *  .  8AT  *1 

=  B  i^r1  +  —  J  +  kl 


8z 


or 


8x 


a2 
9  P, 


ax 


2  2 
a  p  a  p 

'2  .2  2 

9  p  9  p  a  p 

►  _  ft  9AT 

e  ax 

—r +  — r- +  y ' 

ay  az 

9x2  8x8y  8x8 z 

2  2 
a  p  a  p 

+  J  +  y  \ 

2  2  2 
a  p  a  p  9  p 

x  ^  y  _l  ^z 

►*  _  ft  9at 
■  B  9F 

ay  az2 

9y9x  '2  3y3z 

9y  J 

.2  .2  .2 

a  p  a  p  a  p 

-"  '2"  — t~  — t~  y 

ax  3y  3z 


~2  2  a2 
a  p  a  p  9  p„ 

X  +  _ y  +  _  z 

3z3x  3z3y  .  2 


3z 


=  e 


3AT 

3z 


Now,  if  for  thin  plates 
p  -  ctzAT(x,y) 


then 


9P, 

az 


=  ctAT 


axaz 


3AT 

ax 


a2 

8  pz 
3y3z 


=  a 


3AT 

9y 


,  (HD 


(112) 

(113) 

(114) 

(115) 

(116) 
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Assume  p  and  p  are  not  functions  of  z.  Then 


2  2 
9  p  9  p 


3x 


3y 


.2 
a  P, 


3x 


8  Pv  9AT 

- —  +  a  t— 

9x9y  9x 


9AT 

9x 


(117) 


.2 

9  P„ 


9x 


32P 


L  * 


9y 


x  *%■  x  9 AT 

^7^7“  ^ 


=  6 


9AT 

ST 


(118) 


Simplifying  terms, 


32P, 


92p, 


.2 

3 


(!  +  V)  -r  *  TT  *  v  nsf  ■  (fi  -  1,0  >  5F- 


9x 


9y 


(119) 


92p  92p 

(i  +  y)  — 7  +  — T~  + 

9y  9x 


92P  ' 

US?  =  (e  -  ^  3y 


9AT 


(120) 


J.  Thermal  Stress  Deformation  of  a  Thin  Circular  Plate  With  a 
Modified  Correction  for  the  z  Component  of  Deformation 


The  problem  associated  with  assuming  that  the  pz  component  of 
deformation  obeys  Eq.  (113)  is  that  the  p^  contribution  of  deformation 
to  the  p^  component  is  not  accounted  for.  The  following  derivation 
assumes  that  the  material  comprising  the  plate  is  largely  incompressible 
and  that  changes  in  p^,  may  be  directly  related  to  changes  in  the  p^ 
component.  The  coordinates  for  a  circular  plate  are  shown  in  Figure  9. 
For  the  thin  circular  plate,  at  r  the  deformation  is  p^  and  at  r  +  Ar 
the  deformation  is 


(121) 
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Figure  9.  Coordinates  for  a  circular  plate. 
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The  initial  volume  of  the  material  contained  between  r  and  r  +  ir  is 


=  [il  (r  +  Ar)2  -  nr2]  t  , 


(122) 


where  t  is  the  plate  thickness.  Upon  heating  or  cooling  the  plate,  the 
new  volume  is 


r  1 

(  3Pr  1 

I2  „  /  \2" 

vf  = 

*1 

(r  +  Ar  +  Pr  +  57-  Arj 

f  -  "  Vr  *  pr) 

(123) 


where  t  is  the  new  plate  thickness. 
Equating  terms, 


II  ^r  +  pr*  +  Ar^  -  II  (r  +  p^  j  t  =  jn  (r  +  Ar)  -  Ilr 


9pr 

p  '  =  p  +  tt—  Ar 

Kr  r  3r 


p  p  2  2l' 

+  p  '  +  Ar  +  2rp  '  +  2rAr  +  2p  'Ar  -  r  -  p  -  2rp  r  t 

r  r  r  r  r ) 


lr2  +  2rAr  +  Ar2  -  r2  f  t 


V 124  ) 


Let 


P 


r 

r 


p  +  e 
r 


where  e 


(125) 
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Then 


,2 


(pr  *  «) 


2  _  2  2  0 

p  +  e  +  2ep 

r  r 


(126) 


Substituting  terras, 


(e 


t  ^e2  +  2epr  +  Ar2  +  2re  +  2rAr  +  2prAr  +  2Are^ 


^2rAr  +  Ar2} 


(127) 


t 

T“ 


2rAr  +  Ar 


8r 


■Ar 


+  21aF"Ar}pr  +  Ar2  +  2r 


^pr  1  2  ^Pr 

3r  j  r  3r 

(128) 


Simplifying  and  neglecting  small  terms, 


I 


8Pj 

5r  ar 


+  r 


3p 

i 

3r" 


+  r  +  p 


1  + 


3p 

i 

ar" 


Now, 


! 


And  at  z  =  t/2,  the  deformation  is  t  /2  -  t/2. 


(129) 


(130) 
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Using  the  correction  term  of  Eqs.  (114)  and  (130), 


Pz  (r, z)  =  azAT(r)  + 


1  +' 


(9P, 

3r 


-  z 


(131) 


3pT 


Pz  (r,z) 


Z  -  Z  ~  2 


azAT(r)  +  - 5 - 


1  + 


3r 


z[aAT(r)-^l  -(132) 


The  thermostress  equations  in  cylindrical  coordinates  are 


(l  +  y) 


92P, 


1  3  p  p 

r  1  r  r 

72“  +  —  —  ‘  “7 

3r  r 


2 

9  Pz  3AT 

Y  3r3z  K  3r 


(1  +  y) 


2 

3  p 


3r 


„£  +  .L 

7T  3r 


+ 


92P, 

y  =  b 


3AT 


3r3z  3r 


(133) 


Substituting  Eq.  (132)  into  Eq.  (133), 


(1  +  y) 


2  ,  ' 

1  ' 

>+  9  - 

3p 

aAT _ £ 

>  .  B  2*1 

TF  8r  \r  / 

Y  3r 

afli  3r 

k  i 

3r 

(1  +  y)  1 


~2  1 
3  p 

r 

►  4.  < 

(Ml 

3r2 

T 

3r 

t  Y  '  i 

(l  +  Y)  - 


92p„ 


Y  -rr +  Ya  3r~ 

3r 


8 AT  8AT 


8r 


ir + (1  +  y)  h  {r} 


=  (3  -  ay) 


3AT 

3r 


(6  -  ay) 


3AT 

3r 


3_ 

3r 


3p  p 

3^(1  +  .)/ 


(134) 
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Integrating, 


AT 


B 


1 _ 

-  ay 


(135) 


Suppose  the  plate  is  infinite  in  extent;  i.e.,  §  r  =  00 ,  p  -  0, 
Dp^/Sr  =  0,  and  AT  =  0.  This  implies  that  C  =  0.  Therefore, 


AT 


(136) 


for  an  infinite  plate. 

K.  Least  Squares  Method  of  Differentiating  Experimental  Data 

Assume  that  N  measurements  of  radial  deformation  were 

taken  in  some  region  of  a  thin  circular  plate.  Assume  that  the  deforma¬ 
tion  in  this  region  is  given  by 


=  At  +  Br  +  Cr  +  D 


(137) 


Forming  the  square  of  the  difference  (difference  function)  between  the 
curve  given  by  Eq.  (137)  and  the  experimental  data  [5], 


-k  M 


Ar^' 


Br. 

i 


+  Cr.  +  D 


(138) 


Differentiating  with  respect  to  the  unknows  to  obtain  an  extrema, 
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36.  Nr  ,  p  "I 

=  0  =  E  p.  -  At.  -  Br.  -  Cr,  -  D  r  ' 

3A  |_  1  i  i  1  J  1 

36.  Nr  o  ?  "I  ' 

^  =  °  =  i?i  L1  "  AV±  ”  Bl>i  ’  Cri  "  DJ  ^ 

36.  NT  o  o  1 

3CT  =  °  =  .?!  Lpi  -  -  Bri  -  Cri  ’  Dj  ri 


36.  N 

r  -  0  -  I 

i=l 


SD 


-  Ar . 3  r-  Br . 2  -  Cr.  -  D 
i  x  l  i  J 


(139) 


N 

Let  E  =  E 
1=1 


and  rearrange  terms  to  obtain 


A2r.6  +  BSr^  +  CSr..4  +  DSr.j3  =  Sp^r..3 


A2r.^  +  B2r.4  +  C2r.3  +  D2r.2  =  2p.r.2 

1  1  1  1  11 


A2r.4  +  B2r.3  +C2r.2  +  D2r.  =  Sp.r. 
1  1  1  1  11 


A2r.3  +  B2r.2  +  C2r.  +  DN  =  Sp . 

ill  i 


(140) 
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Let  R6 

T  ^ 

=  Er. 

S3 

=  EP^J 

R5 

=  Eri5 

S2 

=  Epjl-j 

R4 

*  *>i4 

SI 

■  Epiri 

R3 

'  i 

S0 

=  EPj 

R2 

v  2 

=  Er . 

R1 

=  Ir. 

R0 

=  N 

Then  S3  =  A*R6  +  B*R5  +  C*R4  +  D*R3 
S2  =  A*R5  +  B*R4  +  C-R3  +  D-R2 


(141) 


SI  =  A*R4  +  B*R3  +  C*R2  +  D*R1 

S0  -  A-R3  +  B*R2  +  C-Rl  +  D.R0  .  (142) 


Equation  (142)  may  be  expressed  in  matrix  form  as 


R6 

R5 

R4 

R3 

~A 

S3 

R5 

R4 

R3 

R2 

0 

B 

S2 

R4 

R3 

R2 

R1 

C 

SI 

R3 

R2 

R1 

R0 

D 

S0 

(143) 


Equation  (143)  may  be  solved  for  A,  B,  C,  and  D— thus  giving  the 
equation  for  with  a  least  squares  curve  fit.  These  curves  can  then 
be  more  readily  differentiated  in  the  thermostress  equations  to  obtain 
better  results , 
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L.  Polynomial  Approximation  to  the  Thin  Plate  Thermostress 
Equation  in  Cylindrical  Coordinates 


The  thin  plate  thermostress  equation  in  cylindrical 
coordinates  may  he  given  as 


AT(r) 


r 


2 

a 


(144) 


The  difference  in  temperature 
given  as 


T  ^  between  location  r  and  R  is 
r-R  o 

o 


(145) 


Suppose 

p  -  Ar2  +  Br  +  C  U46) 

Hr 

for  a  thin  plate.  Then 

=  2Ar  *  B  ,  (H7) 

3r 

and 
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AT 


r-R 


— f  j  2 At  +  B  +  At  +  B  +  —  J 


R 


AT 


r-R 


-T  \3kr  +2B  +  ^ 


R 


ATt-R  =  "r  3ARo  +  2B  +  g-  -  3Ar  -  2B  -  | 
o  l  o 


(148) 


But  C  =  0,  since  p  =  0  at  r  =  0,  Therefore 
r 


4Tr-EQ  *  -3Ar  [Ro  -  r] 


AT 


r-R 


-  3Ar  [r  -  Ro] 


(149) 


Now, 


Ml 

3r 


r  pc.  *  h: 
9r  hr  r 


(150) 


Suppose 


Ar  +  Br  , 


(151) 


then 


r  1“  (2Ar 

9r  l 


+  B  +  Ar  +  B 


)  - 


8  AT 
8r 


(152) 


or 
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9AT 

9r 


(153) 


=  3AT 

Equation  (153)  may  be  used  to  estimate  thermal  gradients  from 
experimental  data. 

III.  EXPERIMENTAL  EXAMPLES 

The  following  three  cases  were  used  to  test  the  theory  presented 
in  Section  II.  Two  problems  are  treated  in  this  section: 

a.  The  heating  of  a  thin  circular  flat  plate  at  its  center  by 
external  heat  generation. 

b.  The  uniform  heating  of  a  rod. 

The  thin  circular  flat  plate  problem  is  treated  in  Sections  IIIB  and 
HID,  and  the  uniform  heating  of  a  rod  is  treated  in  Section  IIIC. 

A.  Computer  System  for  Analyzing  Laser  Speckle  Interferograms 

Figure  10  illustrates  the  automated  system  used  to  analyze 
laser  speckle  interferograms.  A  description  of  the  system  and  its 
functional  components  are  documented  in  reference  [31-  The  system  was 
used  in  the  manual  mode,  and  the  computer  programs  used  to  control  the 
system  are  listed  in  the  appendix. 

B.  Temperature  Change  Measurements  for  a  Heated  Circular  Flat 

Plate  ( Case  I ) 

This  test  was  conducted  to  examine  the  accuracy  in  measuring 
a  temperature  change  AT  in  a  circular  flat  plate  using  laser  speckle 
interferometry.  The  plate  may  be  classed  as  a  thin  circular  flat  plate 
subject  to  a  heat  source  located  at  the  center.  The  heat  source  was 
obtained  from  electrical  resistive  heating  located  at  the  center  of  the 
plate.  The  resulting  temperature  profile  of  the  plate  tends  to  be  a 
function  of  the  radial  distance  from  the  center  of  the  plate. 

Figure  11  illustrates  the  basic  laboratory  configuration  for 
measuring  temperature  changes  in  circular  flat  plates  subject  to  a 
radial  symmetric  temperature  distribution.  In  the  test,  an  18.0-inch 
diameter,  0.125-inch-thick  aluminum  plate  ( Figure  12)  was  mounted  onto 
a  heater  element  support  (Figure  13).  The  heater  element  support  was 
also  aluminum  and  was  covered  with  a  0.25-inch-thick  layer  of 
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Figure  10.  System  used  to  analyze  laser  speckle  interf erograms . 
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Figure  12.  Circular  flat  plate. 
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Figure  13.  Heater  element  support. 
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compressed  K-wool.  Nichrome  wire  was  wrapped  around  the  circumference 
and  length  of  the  support  to  serve  as  a  heater  element.  Again,  K-wool 
was  wrapped  around  the  wire  to  electrically  insulate  it  from  any  metal. 
The  plate  and  heater  were  then  inserted  into  the  plate  heater  unit 
holder  (Figure  14)*  In  order  to  provide  electric  power  access, 
0.13-inch-diameter  holes  were  provided  in  the  holder.  The  plate  heater 
holder  was  then  bolted  to  the  heater  unit  support  base  (Figure  15), 
which  was  then  bolted  to  an  NRC  air  table  for  laser  speckle  work. 

In  order  to  determine  the  accuracy  in  measuring  temperature  fields 
with  laser  speckle,  a  series  of  .0625-inch-diameter,  .095-inch-deep 
holes  were  drilled  at  various  radial  distances  on  the  circular  plate 
(Figure  12).  Copper  ^-constant  on  thermocouples  were  attached  in  each  of 
these  holes  using  a  conductive  epoxy.  This  compound  is  available  from 
Ablestik  Laboratories,  833  West  l82d  Street,  Gardena,  CA  90248,  and 
consisted  of  the  following: 

.  Ablebond  163-4,  part  C  pure. copper  powder 
( 4  parts ) 

.  Ablebond  163-4,  part  A  resin  component 
( 2  parts ) 

.  Ablebond  163-4,  part  B  hardener  component 
(1  part) 

Figure  16  illustrates  a  typical  thermocouple  attached  to  the  aluminum 
plate.  It  is  unfortunate  that  a  significant  amount  of  heat  can  be 
conducted  away  from  the  thermocouple  attachment  point  with  this  design, 
but  it  was  the  best  known  alternative  available.  A  better  alternative 
would  have  been  to  use  a  semiconductor  thermal  radiation  detector  or 
liquid  crystals,  but  these  were  not  readily  available.  Figure  17 
illustrates  the  complete  assembly  of  eight  thermocouples  which  were 
mounted  on  the  circular  flat  plate.  Each  of  the  thermocouples,  includ¬ 
ing  a  reference  at  0°  C,  was  individually  switched  to  a  Leeds  and 
Northrup  8690  potentiometer,  using  the  electrical  skematic  shown  in 
Figure  18.  The  rotary  switch  is  shown  in  Figure  19.  A  voltage  divider, 
using  1-kn  resistors,  was  used  to  plot  thermocouple  output  versus 
thermocouple  number.  The  reference  thermocouple  tempera ture  at  0°  C 
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Figure  14.  Plate  heater  unit  holder 
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Figure  15.  Heater  unit  support  base. 


Figure  16.  Attachment  of  iron-constantan  thermocouples  to 
aluminum  plate. 
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Figure  17. 


Complete  assembly  of  thermocouples  on  the 
circular  flat  plate. 
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Figure  18.  Thermocouple  switch  network  schematic. 
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Figure  19.  Thermocouple  rotary  switch. 
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was  obtained  from  a  reference  bath  ( Figure  20 ) ,  This  device  could 
provide  a  cold  junction  for  several  hours  before  any  ice  replacement 
was  needed.  Figure  21  is  a  photograph  of  the  Leeds  and  Northrup 
Millivolt  Potentiometer. 

In  order  to  supply  electrical  power  to  the  circular  flat  plate 
heater  element  (Figure  11 ),  a  Redlake  Labs  Autotrans former  at  about 
50-V-a.c.  output  (variable)  was  used.  A  piece  of  Nichrome  wire  was 
placed  in  series  with  the  heater  element  to  observe  the  plate  heater 
performance.  When  the  "glow  wire"  was  a  bright  red,  the  proper  current 
flow  was  presented  in  the  heater  assembly.  At  this  point,  the  maximum 
rise  in  plate  temperature  with  time  was  observed. 

Laser  illumination  of  the  plate  was  accomplished  by  expanding  the 
output  of  a  6328  A  spectra  physics  Model  125,  50-mW  He-Ne  laser.  A 
spectra  physics  Model  332  spatial  filter  was  used  for  this  purpose. 
Exposure  control  of  the  laser  was  performed  using  a  Jodon  ES-100 
electronic  shutter  (Figure  22),  The  output  of  the  laser  was  passed 
through  a  horizontally  oriented  cylindrical  lens  to  produce  a  narrow 
concentrated  rectangular  beam  of  light,  which  measured  10.2  mm  by 
228.6  mm  in  the  horizontal  0-direction  on  the  circular  plate.  Beam 
control  was  performed  using  two  X/4  front  surface  mirrors  and  an 
aperture  for  beam  expansion. 

In  the  test,  an  attempt  to  measure  the  temperature  rise  at  station 
#1  (r  =  1.1  in.)  was  made.  Two  32-second  exposures  were  made  using 
AGFA  GAVAERT  holographic  film.  The  initial  station  #1  temperature  was 
81°  F  when  the  first  -  exposure  was  made.  The  plate  was  then  heated  to  a 
thermocouple  measured  temperature  of  112°  F ,  and  a  second  laser 
exposure  was  made.  The  net  temperature  change  at  station  #1  was  31°  F, 
as  measured  with  the  thermocouple. 

Test  Results  -  Case  I: 


S  = 


y-gg  =  2.3195  for  a  9-inch-radius  plate 
74  in. 


81°  F 
112°  F 


measured  with  thermocouple 
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Figure  20.  Reference  junction  ice  bath. 
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Figure  21.  Leeds  and  Northrup  8690  Millivolt  Potentiometer. 
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Figure  22. 


Jodon  ES-100  electronic  shutter  system. 
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Table  1  illustrates  the  data  obtained  from  the  interferogram.  For 
this  data,  r'  is  the  radial  distance  on  the  interferogram  and  r  is  the 
radial  distance  on  the  plate.  x'(r)  is  the  measured  fringe  spacing. 

The  sign  associated  with  x'(r)  determines  the  direction  of  plate  trans¬ 
lation  upon  heating;  (+)  refers  to  motion  in  the  positive  X  direction. 
Table  2  shows  the  data  of  Table  1  corrected  for  true  relative 
translation,  with  the  center  of  the  plate  having  zero  translation. 

Figure  23  illustrates  the  variation  of  y(r)  with  r.  In  order  to 
differentiate  y(r)  with  respect  to  r,  a  quadratic  curve  fit  of  Figure  23 
is  utilized. 

Y(r)  =  =  Ar  +  Br2  +  C  .  (154) 

The  data  for  Eq.  (72)  may  be  tabulated  as 


r 

y(r ) 

.30 

0 

1.23 

.07662 

2.44 

. 19012 

Solving  the  simultaneous  linear  equations  in  A,  B,  and  C  results  in 


A  =  .07415 

B  =  .0054 

C  =  -.02273 

and 

y(r)  =  .07415r  +  .0054r2  -  .02273 
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TABLE  1.  CASE  I  INTERFEROMETRIC  DATA 


r 

X'(r) 

1/x ’ ( r ) 

.13 

.30 

-2.86 

-.34965 

1.05 

2.44 

-3.68 

-.27173 

.53 

1.23 

-3.48 

-.28735 

1,63 

3.78 

-4.68 

-.21367 

2.19 

5.08 

-7.02 

-.14245 

3.44 

7.98 

+  5.75 

+.17391 

3,88 

9.00 

+3.93 

+.25445 

TABLE  2.  CORRECTED  CASE  I  INTERFEROMETRIC  DATA 


r 

1/X(  r ) 

y(r) 

.30 

0 

0 

2.44 

,07792 

. 19012 

1.23 

.06230 

.07662 

3.78 

. 13598 

. 51400 

5.08 

.20720 

1.05257 

7.98 

.52356 

4,17800 

9.00 

.60410 

5.43690 
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Figure  23.  Case  Iy(r)  versus  r. 
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Now  from  Eq.  (55), 


B 


XfS  jl  -  y 
2  a  (  y 


where  A 

f 
S 

a  (Aluminum) 
y  (Aluminum) 


63281  =  2.4913  x  10-5  in. 
74.0  in. 

2.3195 

24  x  10-6/°C 

.3 


Now, 


3 


y(r) 

3r 


.07415  +  . 0108r 


3  =  207.86 


1 

r 


3 


Y(r) 

3r 


|r=l.l  in. 


.0782 


And  from  Eq.  (58), 

AT( r )  =  207.86  (.0782)  °C  s  16.25  °C 

AT(r)  =  29.25  °F  measured  with  laser  speckle  interferometry 
AT(r)  =  31  °F  measured  with  a  thermocouple 


Figures  24  through  27  illustrate  the  laboratory  configuration  without 
the  cylindrical  lens  and  front  surface  mirror  optics  present. 
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Figure  24.  Illumination  direction  for  the  flat  circular  plate. 
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Figure  25.  Side  view  of  Case  I  experiment,  illustrating  the  laser, 
heating  system,  and  thermocouple  read-out  system. 
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0 


Figure  26. 


Case  I  optical  configuration  without  cylindrical 
lens  assembly. 
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Figure  27.  Back  side  of  circular  flat  plate  illustrating  thermocouples. 


77 


C.  Temperature  Change  Measurements  for  a  Heated  Copper  Bus  Bar 
(Case  II) 

In  this  experiment,  a  copper  bus  bar  was  submerged  in  a  heated 
water  bath  and  a  laser  speckle  photograph  was  taken  of  the  heated  rod. 

The  water  was  then  replaced  with  cool  water,  and  a  second  laser  speckle 
photograph  was  taken.  The  double  exposure  interferogram  was  then 
analyzed,  and  the  rod  temperature  was  predicted  based  on  various 
theories  presented  in  Section  II.  The  rod  was  in  a  state  of  uniform 
thermal  contraction  when  cooled. 

The  specimen  was  illuminated  by  a  Spectra-Physics  166  Argon  laser 
operating  at  .9-watt  power.  For  the  test,  the  following  data  were 
observed: 

«  Initial  water  bath  temperature  =  1.280  mV  £.  89.7°  F. 

•  Final  water  bath  temperature  =  .83  mV  =  69.9°  F. 

|  Length  =  9.95  in. 

•  Bus  bar  dimensions  I  Width  =  1.00  in. 

(  Depth  =  0.25  in. 

•  Approximately  5  in.  of  the  bus  bar  was  exposed  in  the  water  bath. 

•  AGFA  10E56  holographic  emulsion  film  was  used  with  two  16-sec 
exposures . 

•  Xf  =  2.4913  x  10-5  (130)  =  .00323869  in.2  (reconstruction  of 

interf erograms  was  with  a  He-Ne  laser ) . 

The  copper  bus  bar  is  a  glass  beaker  immersion  tank,  the  Case  II 
laboratory  setup,  and  a  side  view  of  the  Case  II  experiment  are  shown  in 
Figures  28,  29,  and  30,  respectively. 

The  data  obtained  from  the  interferogram  is  illustrated  in  Table  3. 
For  this  table, 

L  =  distance  from  an  arbitrary  reference  point  on  the 
interferogram 

p  =  longitudinal  specimen  deformation  at  location  L. 

X  =  fringe  spacing  in  the  interferogram  analyzer  plane. 

A £,  =  net  longitudinal  deformation  measured  from  the  arbitrary 

reference. 

All  data  was  measured  in  the  symmetry  plane  of  the  specimen. 
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*  5  "  H-t 


Figure  28.  Copper  bus  bar  in  glass  beaker  immersion  tank;. 

( Note  thermocouple  leads  in  the  tank  used  to 
monitor  water  temperature.) 
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4^:#%  *  * 


Figure  29.  Case  II  laboratory  setup. 

(Note  water  drainage  sys¬ 
tem,  thermocouple  tempera¬ 
ture  read-out  equipment , 
and  the  optical  system. ) 
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Figure  30.  Side  view  of  Case  II  experiment. 
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TABLE  3.  CASE  II  LABORATORY  DATA 


L 

( in. ) 

X 

(  in.  ) 

P 

( in. ) 

Afc  =  p  -  p 

0 

5.54 

.0005846 

0 

.5 

5.00 

.0006477 

6.310  x  10*  5 

1.0 

4.30 

.0007531 

1.685  x  10'4 

1.5 

3.94 

.0008220 

2.374  x  10"4 

2.0 

3.76 

.0008613 

2.767  x  10*4 

For  this  particular  case,  the  temperature  change  can  he  approximated 
from 


AT 


1  A£ 

o  rr 


(155) 


The  results  of  this  case  are  shown  in  Table  4>  where 


TABLE  4.  CASE  II  TEMPERATURE  MEASUREMENT  RESULTS 


L 

( in. ) 

n 

( in.  ) 

A£ 

L 

AT 
(  °F ) 

0 

0 

- 

~ 

.5 

6.310  x  10"5 

1.262  x  10-4 

16.12 

1.0 

1.685  x  10~4 

1.685  x  10"4 

21.52 

1.5 

2.374  x  10~4 

1.582  x  Kf4 

20.21 

2.0 

2.767  x  10"4 

1.383  x  10'4 

17.66 

The  average  temperature  change  computed  using  Eq.  (155)  was 
18.87°  F,  which  is  about  a  4.8  percent  difference  from  the  thermocouple 
reading  of  19.8°  F. 

D.  Temperature  Change  Measurements  for  a  Heated  Circular  Flat 
Plate  (Case  III) 

This  test  was  identical  to  Case  I.  The  data  for  this  test 
case  is  illustrated  in  Table  5. 

TABLE  5.  CASE  I  TEST  DATA 


r  ( in. ) 

Pr  (in. ) 

.30* 

-.001495 

1.23* 

-.001228 

2.44 

-.001161 
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For  the  data  in  Table  5  marked  with  an  asterisk,  assume  that 


p  =  Ar  +  B 
Hr 

Therefore, 

-.001495  =  ,3A  +  B 

-.001228  =  1.23A  +  B  .  ( 156) 

From  Eq.  (156), 

.A  =  2,8709  x  10-4 

B  =  -1,5811  x  10' 3 

B  is  the  offset  error  in  the  data  of  Table  5,  due  to  a  uniform  trans¬ 
lation  of  the  test  fixture.  Therefore,  the  p  component  of  deformation 

-3 

in  Table  5  should  be  corrected  by  1.5811  x  10  in.  to  account  for  this 
uniform  translation.  The  corrected  data,  Pr%  Is  illustrated  in 
Table  6. 


TABLE  6.  CORRECTED  TABLE  5  DATA 


r 

P 

P.J 

r 

r 

.30 

-.001495 

.0000861 

1.23 

-.001228 

.0003531 

2.44 

-.001161 

.0004201 

84 


L 


Now  at  location  r  =  1.23  in., 


F"  =  ~'~"l'  2310 —  =  2.8707  x  10~4 


9p 

1 

8F 


=  3.5J1  x  10^8^61  x  IQ'5  .  2. 8709  x  10~4 


9p 

I 

9r" 


_  4.201  x  10-4  -  3.531  x  10“4 

£.44  -  1.23 


5.5371  x  10 


-5 


9p 

r 

w 


=  1.7123  x  10 


-4 


avg 


From  Eq.  (108),  with  C  =  0  for  an  infinite  plate  solution. 


j3pr  pr 
AT  =  r  {tt^-  +  — 
1 9r  r 


And  from  the  experimental  data. 


I1-' 


AT  =  65625  {1.7123  x  HT4  +  2.8707  x  10~4!  °F 


AT  -  30.075°  F  at  r  =  1.23  in. 


AT  ~  31°  F  at  r  =  1,10  in.  (thermocouple  measurement) 


From  the  above  results,  agreement  is  very  good. 


(157) 
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IV.  CONCLUSIONS 


The  theory  for  laser  speckle  interferometric  thermoelasticity  was 
presented  in  Section  II.  Two  specific  problems,  and  three  methods  of 
data  analysis  were  presented  in  Section  III.  Laser  speckle  interferom¬ 
etry  may  he  used  to  measure  the  in-plane  deformation  of  a  body,  hut  no 
information  about  the  out-of-plane  component  is  obtainable.  This 
necessarily  implies  that  certain  simplifications  to  the  theory  are 
needed.  After  making  estimations  for  the  out-of-plane  component  of 
deformation  for  a  body,  the  general  thermoelasticity  equations  are 
solvable , 

Section  III  illustrates  that  for  the  case  of  thin  plates  and 
uniformly  heated  rods,  accurate  temperature  measurements  can  he  made. 

As  long  as  the  deformation  of  some  region  in  a  body  is  attributed  to 
localized  heating,  then  accurate  temperature  measurements  can  he  made. 
Major  problems  occur  when  the  deformation  is  attributed  to  thermostress 
generated  elsewhere  in  a  body.  This  case  occurs  at  the  outer  radius  of 
the  circular  flat  plate  where  the  deformation  and  derivatives  of 
deformation  are  not  zero.  This  is  the  result  of  thermostress  generated 
at  the  center  of  the  plate  and  is  very  difficult  to  analyze. 

The  thermostress  equations  usually  require  some  form  of 
simplification,  for  they  involve  partial  derivatives  of  the  deformation 
field.  After  a  suitable  set  of  approximations  are  made,  the  thermo¬ 
elasticity  equations  may  he  made  more  amenable  to  analysis.  Thermal 
gradient  measurements  are  easily  made,  since  many  approximations  to  the 
deformation  and  temperature  fields  are  not  required. 
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APPENDIX 


Computer  Codes  Used  to  Analyze 
Laser  Speckle  Interferograms 
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Computer  Code  -  1 


C - CIRCULAR  LASER  SPECKLE  INTERFEROMETRY  READER 

C - WRITTEN  BY:  JOHN  A  *  SCHAEFFEL f  !R. 

DIMENSION  XY ( 10 » 50 ) »XX < 10 »S0 ) * YY< 10 » 50 ) 

WRITE ( 5  f 1 ) 

1  FORMAT  <  '  INPUT  R1  rRNrTl  »TM»NfM?  IS--4F5.0»3I2'  ) 

READ ( 5  > 2 )  R1 » RN » T1 » TM » N  »M r IS 

2  FORMAT (4F5.0»3I2) 

WRITE ( 5  j  3  ) 

3  FORMAT ( '  INPUT  SF , E-2F10 . 0 ' ) 

READ ( 5  j  4 )  SF  t E 

4  FORMAT (2F10.0) 

XL=0  ♦ 

YL=0  ♦ 

DR=(RN~R1 > /FLOAT (N~l ) 

DT=<TM-T1 > /FLOAT <M-1> 

DO  7  I=1»M»1 
DO  7  J=1 »N»1 
T=T1+FL0AT  < 1-1 )*DT 
R=R1+FL0AT ( J~1 >*DR 
TP=3.14159*T/180, 

X=  ( R/SF ) *COS ( TP ) 

Y=  ( R/SF ) #SIN ( TP ) 

IXX=INT((X-XL)/.001) 

IYY=INT( <Y-YL>/.001  ) 

IF ( IXX ♦ NE . 0 )  CALL  XADV<IXX»IS) 

IF ( I YY ♦ NE ♦ 0 )  CALL  YADV<IYY,IS> 

XL=XL+FLOAT ( IXX)/1000* 

YL=YL+FLOAT ( I YY ) /1000 ♦ 

WRITE(5»5) 

5  FORMAT ( '  DATA-2F5 . 0  J ' ) 

READ <5, 6)  XS  t FS 

6  FORMAT (2F5.0) 

XY(I» J)=SF*E/XS 

XX(I»J)=XY(I » J) *COS  ( FS*3 . 14159/180  .  ) 

YY ( I » J )-XY ( I » J )*SIN (FS*3 .14159/180, ) 

7  CONTINUE 

DO  9  1  =  1 » M » 1 
DO  9  J=1 » N » 1 
T=T1+FL0AT ( I~1 >*DT 
R= ( Rl+FLOAT ( J-l )*DR) 

WRITE ( 5 » 8)  I» J»T»R»XY(I» J) »XX< I » J) »YY(I r J) 

8  FORMAT ( '  I~\r  I3»3X* ' J=' r I3»3X» ' ANGLE® ' ? F8 . 2* 3X t 'R= ' »F8 . 2 1 3X t 
1  'XY='  ,F10.6»3X»  'XX= '  »F10  >  6»3X»  ,YY-~/F10.6) 

9  CONTINUE 
IXX=~INT ( XL/ ♦ 001 ) 

I YY=-INT ( YL/ ♦ 001 ) 

IF ( IXX ♦ NE « 0 )  CALL  XABV< IXX , 10 > 

IF ( I YY  *  NE  « 0 )  CALL  YADUC TYY , 10 > 

STOP 

END 
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rj  CJ  f\  tH  C-4  n  00  ■‘x"  iD  -O  O  O 


Computer  Code  -  1  (Continued) 

S U B R 0 U T I N E  Y A D Vi  IS ? I ft > 

- IS “NO *  STEPS  ( T-FWB y -=REV ) 

- IR=ADVANCE  RATE  OF  STAGE 

x=o* 

IF ( IS ♦ GT ♦  0 )  GOTO  3 
IF-IABS(IS) 

DO  2  I  “  1  y  I P  y  1 

.  CALL  IP0KE< “167772 y “020000) 
DO  7  K— 1 y IR  y 1 
Y-SIN(X) 

CALL  I POKE ( *  167772  y  *  000000 ) 
DO  1  J=1 y IR  y 1 
Y-SIN ( X ) 

CONTINUE 

GOTO  6 

CONTINUE 

DO  5  1 1  —  1 y I S  y 1 

CALL  IP0KE( “167772 y  *010000) 

DO  8  KK-1 y IR y 1 

Y-SIN ( X ) 

CALL  IP0KE( n167772y "000000) 
DO  4  JJ-lylRyl 
Y-SIN ( X ) 

CONTINUE 

CONTINUE 

RETURN 

END 

SUBROUTINE  XADV(ISylR) 

- I S— NO ♦  STEPS  C+=FWDy-=REO) 

- IR-ADVANCE  RATE  OF  STAGE 

X=0* 

IF( IS ♦ GT ♦ 0 )  GOTO  3 
IP=IABS ( IS ) 

DO  2  I  —  1 y I P  y 1 

CALL  I POKE ( “1 67772 y “100000) 
DO  7  K-l y IR y 1 

7  Y=SIN ( X ) 

CALL  IPOKE( “167772? “000000) 
DO  1  J-l ? IR  ? 1 

1  Y-SIN ( X ) 

2  CONTINUE 
GOTO  6 

3  CONTINUE 

DO  5  II-l yISyl 

CALL  I POKE ( “ 1 67772  y “ 040000 ) 

DO  8  KK=1 ? IR ? 1 

8  Y-SIN ( X ) 

CALL  I POKE ( “167772 ? ” 000000 ) 
DO  4  J J-l y IR  y 1 

4  Y=S IN ( X ) 

5  CONTINUE 
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w  ro  *-*  n  n  n  n 


Computer  Code  -  2 

CYLINDRICAL  COORD INA I  Lb 
LASER  SPECKLE  INTERFEROMETRY 
TEMPERATURE  CALCULATOR  CODE 
WRITTEN  BY  JOHN  A.  SCHAEFFEL  t  JR. 

DIMENSION  A  ( 2  i> 2 ) 

WRITE(5r 1 ) 

FORMAT ( '  INPUT:  SF > E > PSI-3F10 ♦ 0 ' ) 

READ  ( 5  *  2 )  SFrErSI 
FORMAT ( 3F 10.0) 

WRITE ( 5  r  4 ) 

FORMAT ( '  INPUT  R  ?  THETA  * DR-3F10 ♦ 0  '  ) 

READ  (Sr  2 ) 

X=R*C0S<T*3. 14159/180. )-D 
Y=R*SIN<T*3. 14159/180. > 

IX=INT<X*1000./SF> 

IY=INT<Y*1000./SF> 

IR=INT<D*1000./SF> 

IRR=2*IR 

IF(IX.NE.O)  CALL  XADV(IXs-5> 

IF(IY.NE.O)  CALL  YADV< IY»5> 

DO  7  I  — 1  *  2  r 1 
WRITE(5*5) 

FORMAT ( J  INPUT  HORIZONTAL  FRINGE  SPACING-F5 . 0 J ' > 
READ <5 >6)  DH 

6  FORMAT (F5.0) 

A<lrI)=SF*E/DH 
CALL  XADV < IRR  r  5  > 

7  CONTINUE 
IB=-3*IR 

CALL  XADV < IB» 5) 

IRRR=-IR 

CALL  YADV ( IRRR r 5 ) 

DO  10  1  =  1 r  2  *  1 
WRITE ( 5  r  8 ) 

8  FORMAT  (  '  INPUT  VERTICAL  FRINGE  SPACING-F5 . 0  .*  '  ) 

READ  <  5*  9  >  DO  . 

9  FORMAT <  F5 ♦ 0 ) 

A ( 2r I ) =SF*E/DV 
CALL  YADV(IRRr5) 

10  CONTINUE 

CALL  YADV< IB»5) 

IX=-IX-IR 
I Y=-I Y 

IF(IX.NE.O)  CALL  XADV(IXr5) 

IF(IY.NE.O)  CALL  YADV< 1Y »5> 

DT=SI* ( A< 1 1 2 ) -A ( 1 » 1 ) +A  <  2  *  2 > -A  <  2 1 1 ) ) / <  2 . *D ) 

WRITE <5>11)  DT 

11  FORMAT (  '  DELTA  TEMPERATURE  CHANGE== ' »F6. 1 ) 

GOTO  3 

STOP 

END 

SUBROUTINE  YADV<IS»IR> 
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in  -o 


Computer  Code  -  2  (Continued) 

c - IS=NO.  STEF'S  <+=FWD»-=REV> 

c - I R- ADVANCE  RATE  OF  STAGE 

X=0. 

IF< IS*GT»0)  GOTO  3 
IF-IABSCIS) 

DO  2  I=1?IP»1 

CALL  IF'OKE  (  "  167772  ?  “  020000 ) 
DO  7  K=1*IR»1 

7  ’  Y-SIN ( X ) 

CALL  I POKE ( " 167772 » “ 000000 ) 
DO  1  J=l»  IRs-1 

1  Y=SIN ( X ) 

2  CONTINUE 
GOTO  6 

3  CONTINUE 

DO  5  11=1 r IS* 1 

CALL  IF'OKE  (  "  167772  t  n0J  0000) 

DO  8  KK=1»IR?1 

8  y=SIN ( X ) 

CALL  IF'OKE  <  ’  167772  ?  “  000000  ) 
DO  A  J J=1  > 1 1 R  f 1 
A  Y=SIN<X> 

5  CONTINUE 

6  CONTINUE 
RETURN 
END 

SUBROUTINE  XADV<IS>IR) 

C - IS=NO.  STEF’S  <  +=FWD * ~=REV> 

C - IR=ADVANCE  RATE  OF"  STAGE 

x=o  ♦ 

IF ( IS  *  GT ♦ 0 )  GOTO  3 
IF'=IABS<  IS) 

DO  2  I=1»IP»1 

CALL  IF'OKE (  "  167772 » ’  100000) 
DO  7  K=1 » IR » 1 

7  Y=SIN(X) 

CALL  IF'OKE  (  "  167772  t  "  000000  ) 
DO  1  J=lfIR»l 

1  Y=SIN ( X ) 

2  CONTINUE 
GOTO  6 

3  CONTINUE 

DO  5  II=l»ISfl 

CALL  l'POKE  C  “  167772  r  11  0-40000  ) 

DO  8  KK-l » IR  y 1 

8  Y-SIN ( X ) 

CALL  IF'OKE  <  "167772  »  ”000000) 
DO  A  JJ=lfIR»l 
Y=SIN<X) 

CONTINUE 
CONTINUE 
RETURN 
END 
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-fc>  w  r  o 


Computer  Code  -  3 

C - CYLINDRICAL  COORD  I  NAT  t.  S 

C - LASER  SPECKLE  INTERFEROMETRY 

C - TEMPERATURE  CALCULATOR  CODE 

C - WRITTEN  BY  JOHN  A.  SCHAEFFEL  r  JR. 

WRITE (5*1) 

.1  FORMATS  INPUT:  SF,E»PSI-3F10.0'  ) 

READ <5* 2)  SFfEfSI 
FORMAT (3F10.0) 

WRITE (5*4) 

FORMAT ( '  INPUT  R r THETA t DR-3F5 .O') 

READ ( 5  y  5 )  R  >  T  »  D 
FORMAT (3F5.0) 

XI =<R-D)*CQS<T*3. 14159/180. ) 

Y1=(R-D)*SIN<T*3. 14159/180.) 

X2==  ( 2. *D)*CQS<T*3. 14159/180.  ) 

Y2= < 2. *D)*SIN(T*3. 14159/180. > 

IX1=INT(X1*1000./SF> 

I Yl-INT ( Yl*l 000 ♦ /SF ) 

IX2-INT ( X2*1000 ♦ /SF ) 

IY2=INT(Y2*1000./SF> 

IF(IXl.NE.O)  CALL  XADU(IX1»5) 

IF(IYl.NE.O)  CALL  YADU(IYi»5) 

WRITE (5*6)  T 

6  FORMAT ( '  FRINGE  SPACING  AT' r IX rF5 . 1 t IX »' DEGREES ' ) 
READ<5»7)  DH 

7  FORMAT (F5.0) 

A=(R-D)*SF*E/DH 

IF ( 1X2 ♦ NE . 0 )  CALL  XADU<IX2»5> 

IF ( IY2 . NE . 0 )  CALL  YADU<IY2,5> 

WRITE  <  5  y  8 ) 

8  FORMAT ('  FRINGE  SPACING  AT 'y IX y F5 . 1 r IX y' DEGREES ' ) 
READ ( 5  y  7 )  DH 

B=<R+D)*SF*E/DH 

C=( ( <B-A)/<2.*D> )/R)*SI 

WRITE(5y9)  C 

9  FORMAT ("  DELTA  TEMPERATURE  CHANGE^ ' y F5 . 1 ) 

IX=- ( IX1+IX2) 

IY=-( IY1+IY2) 

IF(IX.NE.O)  CALL  XADU<IXy5> 

IF(IY.NE.O)  CALL  YADU < IY  y  5 ) 

GOTO  3 

STOP 

END 
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no  i>  cn 


Computer  Code  -  3  (Continued) 

SUBROUT INE  YADV <  IS » IR ) 

C - 1 S ~N0 .  STEPS  (  +=FWD  r  ~=REU  J 

C - I R= ADVANCE  RATE  OF  STAGE 

x=o  ♦ 

IF(IS.GT.O)  GOTO  3 
IP-IABS ( IS ) 

DO  2  I=lrIP»l 

CALL  I POKE ( " 167772 1 ”020000) 
DO  7  K=1»IR»1 

7  Y=SIN(X) 

CALL  IPOKE ( “ 167772  *  *  000000 ) 
DO  1  Jsl'IR,l 

1  Y-SIN ( X ) 

2  CONTINUE 
GOTO  6 

3  CONTINUE 

DO  5  11=1 ? IS >1 

CALL  IPOKE ( “167772? ”010000) 

DO  8  KK=1 » IR r 1 

8  Y=SIN ( X ) 

CALL  IPOKE < "167772 r ”000000) 
DO  4  JJ=l»IRfl 

4  Y=SIN ( X ) 

CONTINUE 
CONTINUE 
RETURN 
END 

SUBROUTINE  XADV(ISrlR) 

- IS=NO ♦  STEPS  <  +=FWD  ? -=REV ) 

- IR=:ADVANCE  RATE  OF  STAGE 

X=0  ♦ 

IF  ( l'S  ♦  GT .  0  )  GOTO  3 
IP=IABS( IS) 

DO  2  1  =  1 r IP  t 1 

CALL  IPOKE ( " 167772  r ”100000) 
DO  7  K=1 y IRf 1. 

7  Y=SIN  <  X ) 

CALL  IPOKE ( “167772? ”000000) 
DO  1  J=1 ?  IR?  1 

1  Y=S1'N  (  X ) 

2  CONTINUE 
GOTO  6 

3  CONTINUE 

DO  5  11=1 ? IS, 1 

CALL  IPOKE < ” 167772 r ”040000) 

DO  8  KK=1 ? IR ? 1 

8  Y =SIN ( X ) 

CALL  IPOKE < " 1 67772 » ” 000000 ) 
DO  4  JJ=1 ?  IR  ?  1 

4  Y=SIN ( X ) 

5  CONTINUE 

6  CONTINUE 
RETURN 
END 
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LIST  OF  SYMBOLS 


A,B, C  y D 

Constants 

B. 

1 

Body  Force  Per  Unit  Volume 

d,dH'dv 

A  A  A 

e  .e  ,e 
x'  jr  z 

AAA 

e  ,e.,e 
r’  e’  z 

Fringe  Spacing  From  Laser  Speckle  Interf erograms 

Cartesian  System  Unit  Vectors 

Cylindrical  System  Unit  Vectors 

E 

Modulus  of  Elasticity 

f 

Distance  Between  Laser  Speckle  Interferogram  and  the 
Analyzer  Plane 

G 

Lame 1  Constant 

AAA 

Cartesian  Unit  Vectors 

l,m,n 

Cartesian  Direction  Cosines 

L 

Length  of  Specimen 

m 

Fringe  Order 

n  ,n  ,n 
x’  y’  z 

Stress  Plane  Direction  Cosines 

r,  0,z 

Cylindrical  Coordinates 

S 

Film  Scale  Factor 

«A  -A  -A 

Sx,  Sy,  Sz 

Stress  Vectors 

t 

Plate  Thickness 

u,v,w 

Components  of  Deformation  in  Cartesian  Coordinates 

V  UV 

Components  of  Deformation  Determined  by  Laser  Speckle 
Interferometry 

U. 

1 

Components  of  Deformation  in  Cartesian  Coordinates 

x,y,z 

Cartesian  Coordinates 

xj 

Cartesian  Coordinates 

a 

Coefficient  of  Thermal  Expansion 
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LIST  OF  SYMBOLS  (CONCLUDED) 


e 

Y 

r 

Ae 

AT 


e 

e 


ij 

t 


0 

\ 

V 

P 


°n 

♦ 

x(r) 


h 


Q 

— k 

V 


Vo 


V 


2 


Material  Constant 
Material  Constant 
Material  Constant 
Total  Change  of  Length 
Temperature  Change 
Strain  Tensor 
Thermal  Strain 
Angle 

Wavelength  of  Light,  Material  Constant 

Poisson  Ratio 

Deformation  Vector 

Stress  Tensor 

Stress  on  Plane  n 

Strain  Potential 

Fringe  Spacing  at  Radius  (r) 

Material  Constants 
Some  Region  of  a  Body 
Gradi ent 
Divergence 
Laplacian 
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